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Abstract 

We point out that detectable inflationary tensor modes can be generated by particle or string sources 
produced during inflation, consistently with the requirements for inflation and constraints from scalar 
fluctuations. We show via examples that this effect can dominate over the contribution from quantum 
fluctuations of the metric, occuring even when the inflationary potential energy is too low to produce a 
comparable signal. Thus a detection of tensor modes from inflation does not automatically constitute a 
determination of the inflationary Hubble scale. 
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1 Introduction and Motivations 

Gravitational radiation, detectable through B mode polarization in the CMB PQ, provides an important 
handle on primordial cosmology. In the context of inflation, quantum fluctuations of the tensor modes 
in the metric yield a power spectrum of the form 

(hih( f ) = (27rf5(k + k')5 ss 'V h , P h = ~~ , (1) 

whose amplitude is given directly by the scale of the inflaton potential V{<p) ~ H 2 Mp (the index s here 
labels the polarization state) . For this reason, a detection of primordial B modes is often identified with 
a measurement of the Hubble scale H during inflation. It is also related to the range of the inflaton [3] in 
Planck units, providing an ultraviolet-sensitive observable. Observational projects sensitive to B modes 
are expected to test GUT-scale inflation and Planck-scale field ranges in the relatively near term. To be 
detectable at least in near-term observational projects)]] the amplitude must satisfy 

h ~ h k k 3 / 2 > 10- 6 (2) 

In this paper we will consider additional sources of gravitational waves which may be present during 
inflation^] The inflaton might generically be expected to couple to other degrees of freedom X and 
produce excited states of these sectors as it rolls through points in field space where they become light. 
These excited degrees of freedom can source gravitational waves (GWs), which freeze out as they cross the 
horizon. A basic question is to what extent these additional sources can produce detectable primordial 
gravity waves, and how their amplitude hx and other properties compare to those produced by the basic 
process Q. 

Of course the energy density px in the X sector dilutes away during inflation because of the exponen- 
tial expansion. A single production event would lead to a scale-dependent feature in the GW spectrum, 
an interesting possibility in itself. However, during inflation <p may well encounter multiple points with 
new light degrees of freedom. For example, in the mechanism [8] [9] for inflation along angular (axion) 
directions which are typically extended via monodromy, this is automatic: any such production process 
is repeated at short intervals because of an underlying circle in field space. This repetition of the process 
replenishes the supply of extra modes as inflation dilutes them. The requirement of reheating involves 

^ee [2] and references therein for a recent treatment of observational expectations. 

2 See [3] for an interesting analysis of gravitational waves from phase transitions during inflation, and 5 for 
earlier work on some effects of particle production during inflation, in the regime of parametric resonance [5J. 
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coupling of the inflaton to other sectors after exit, and in [5J [9] it would be interesting to study the 
implications of this for earlier particle production events. In other cases with a sufficiently rich spectrum 
of fields, repeated production events may also occur, and in general it is interesting to consider possible 
observational consequences. 

With this motivation, in this note we estimate hx for particles and for strings produced when their 
mass or tension depends sufficiently strongly on the inflaton. We find that the amplitude of the GWs 
generated by various processes including particle/string production, Bremsstrahlung, string oscillations, 
and decays can compete (and in some cases exceed) the contribution Q. This is of interest for two 
reasons: 

• It indicates that even within the context of inflation, the observation of a scale invariant spectrum 
of B modes does not automatically constitute a direct measurement of the inflationary potential. This is 
relevant for the goal of having a systematic treatment of inflationary mechanisms and signatures. 

• It provides a new regime in which to look for observational signatures of exotic sources J^] 

• It introduces an additional (model-dependent) signature of some classes of inflationary mechanisms 
such as [8]0 

The paper is organized as follows. We start in the next section by deriving the basic requirements 
for a detectable signal and how this fits well within the basic bounds on energy densities of extra sectors 
during inflation. Then we review and apply the standard calculation of gravitational waves produced by 
stress energy sources, giving several examples of particle and string sources which generate a competitive 
signal. In an appendix we will review particle production in the presence of a time-dependent mass, 
and describe the salient features we will need of string production in the presence of a time-dependent 
tension. 



2 Basic checks 

A basic requirement of the additional sources is that their energy density px be subdominant to the 
inflationary potential energy: 

px < V ~ H 2 M 2 P (3) 

by at least a factor of the slow roll parameter e ~ jpi ^5 10 _1 - Given that they are produced by the 
rolling scalar field, the energy density px in the additional sources will be at most of the order 

<p 2 ~ eH 2 Mp. (4) 

It will be of this order if the inflaton dumps a significant fraction of its kinetic energy into the X sector, 
as can happen naturally on a steep potential [9j. The fraction / = px/H 2 Mp < e of the total energy 
density which is carried by the sources will figure into our estimates below for the strength of the GW 
signal. 

3 Thc interesting possibility of detecting post-inflationary exotics such as cosmic strings has been studied exten- 
sively; here we are concerned with a distinct window of potential signatures. 

4 This mechanism makes predictions for the amplitude of GWs from tensor fluctuations and for the tilt of the 
spectrum, and depending on the model parameters can lead to additional signatures which are more detailed, such 
as oscillations in the power spectrum and resonant non-Gaussianity [51 lit). The present work introduces another 
model-dependent signature. 
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Once produced by the extra sources X, the additional gravitational waves satisfy the standard equa- 
tion of motion in the inflationary near-de Sitter background, forming linear combinations 



w ^ (• + am) + A * {- + am) '-' tMt)H • (5) 

with A(i 2)k coefficients determined by the initial conditions, which freeze out when k/a(t) reaches the 
Hubble scale H as usual. (This solution ignores the slow time-variation of H.) Gravitational waves that 
are produced inside the Hubble patch decrease in amplitude by a factor of H/(k/a,i) before freezing out, 
where a% is the scale factor at the initial time of production of the mode and hja% ~ UJi its physical 
momentum at that time. 

Gravitational waves of initial frequency OJi make up an energy density at freeze out of order 

PGW ~ h{t) 2 M 2 P \ freezeout ~ ojfM^h 2 (~\ ~ Pi (^) (6) 

where pi is the initial energy density contained in the GWs and hi their amplitude at the time of their 
production; in the last factor we took into account the redshifting noted above of the modes and of the 
frequencies. We are interested in whether the resulting frozen out modes can be competitive with those 
from GUT-scale inflationary theory and visible in near-term observations. Comparison with ^ shows 
that for detectability and consistency with the e condition Q we must require 

lO" 6 < fc- < Vi (-V ^ < 10V/ 4 i? (7) 



OJi \ UJi 



Thus the possibility of observable GWs sourced by the X sector is not immediately excluded by any 
simple consideration of energetics. In order to determine if this possibility is viable, we must work out 
the spectrum of frequencies OJi in concrete examples. In particular, as we will discuss further below, 
low-frequency GWs can be suppressed by interference from multiple scattering events in a dense gas of 
X particles or strings. It is interesting to note that if the X sector degrees of freedom decayed into 

high-frequency GWs with a typical frequency of order \J~<j^ = \/eHMp, the condition (7) translates 

into the condition jg- > 1CP 6 , which is to say that the effect would still be marginally competitive. 

We should also emphasize that the X sector can in general also emit scalar perturbations 5<f); since 
its production arises from its coupling to the rolling inflaton field <f>, at the time it is created there is 
a nontrivial coupling between X and 5<p. Its coupling to 5<j) at later times is model dependent (related 
to the functional form of the ^-dependent mass or tension). In each of the examples below, we will 
determine the strength of the scalar perturbations and estimate their non-Gaussianity to ensure that 
they are consistent with phenomenological constraints. 

It will be useful to rephrase the condition for detectability in terms of the energy density contained 
in gravitational waves. In general, from Q, this requires 

t 2 1 PGW . -in-12 (o\ 

n | freezeout ~ ff2 M 2 - iU W 

and below we will estimate this quantity for particle and string sectors X. 
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3 Gravitational Wave Sources 



In this section, we will start by reviewing the standard derivation of gravitational radiation, following 
the comprehensive treatment in [12], which includes a detailed analysis of Bremsstrahlung. We will wish 
to generalize this analysis in several ways. In particular, we will consider decays and production events 
as well as Bremsstrahlung. Also, we would like to include the effects of the inflaton field 4> coupling to 
particle or string sources, with <j) determining their mass or tension. 



3.1 General setup 

Let us begin by briefly collecting some of the basic results on gravitational radiation, and set up our sys- 
tem. We will shortly make simple estimates in special cases, but it is worthwhile to first lay out the general 
problem. Given the stress-energy of sources, at the linearized level one obtains tensor perturbations [12] 



^ ~ SvrMj, ./ |x - x'| M ' ' + 1 ' 



where S^ u = T^ u — ^T^Tj^x, w). It is convenient to work with the Fourier transform of the stress-energy 
tensor, and one finds a result for the total energy emitted per solid angle 



dE 



dn 8ttMJ, 



J doju 2 (T x »*(k,u;)T Xu (k,oj) - -|T*(k, 



^)l 2 • (io) 



A particularly simple situation to consider is one in which particles of fixed (time- independent) masses 
scatter and emit GWs through Bremsstrahlung. As derived in |12j . that leads to a GW signature with 
total energy 



( dE \ __ ur v - VnVm 
\dncbjj ~ 2vr 2 M2 (P N ■ k)(P M • k) 



(P N ■ P M ) 2 ~ 2 m N m M 



(11) 



where we have taken the limit that the wavelength of the emitted Bremsstrahlung radiation is long 
compared to the scattering time and the time between scattering events. Here N, M index particles with 
momentum Pn,Pm, and ?7jv,?7m are il depending on whether the particle is ingoing or outgoing in a 
given event. In particular, if there is no scattering, so that the incoming and out going momenta are the 
same, (11) gives zero via the cancelations arising from the r]N,ilM factors. 

For processes occurring well within the Hubble scale H~ , stress-energy is conserved to a good 
approximation, shared between particle or string sources and the inflaton field. Before considering specific 
examples, let us briefly set up the full problem. The classical action is 

d A x^[] i M 2 P n + CA+Sx + SxY + SY (12) 



Sx = ~J2 I I rfT^(^-x^(r))m(0(t,x))^-g^(x p (r)) "^^"^' / g(t-t p ) (13) 



v 



d^x j "dr^Cx" -x£(T))m(0(i,x))^ ~ <~ dx " (r) ' 



J d 4 x J d 2 a5^(x» - x^( < 7))r(0(t,x)) v /-Detg^(x s ( < 7))a Q x^( < 7)^x"(r)fl(t -t s ) 
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where p indexes particle sources and s indexes string sources with a mass m or tension T which in general 
depend on the inflaton field eft. The step functions in ( 12 ) reflect the fact that the particle or string sources 
are produced during inflation, on a timescale short compared to the scale of gravitational waves we wish 
to consider. 

Here Sxy describes couplings of the X sector (particles or strings) to other fields Y, included because 
this is a generic possibility which affects their decays. Of course the X sector necessarily couples to gravity, 
and in addition to producing GWs classically, X strings can decay into gravitons, and two or more X 
particles can annihilate into gravitons. 

Clearly, solving for the detailed dynamics and GW spectrum from (12) is prohibitively difficult in 
general since the system is nonlinear; but this is also not necessary for our goal of estimating the leading 
contributions in some cases. For sufficiently large density, for example, one may instead treat the collection 
of sources X as a fluid, along the lines recently reviewed and applied in [13]. For sufficiently small density, 
it is tractable to sum the effects of individual sources. 



3.2 Examples of competitive effects 

Let us now consider some illustrative cases where the GW signature is competitive with (or exceeds) 
the tensor modes arising from the standard mechanism ([!]) . We will consider the production event itself, 
the effect of decays of the produced particles, and the effect of ordinary Bremsstrahlung radiation from 
scattering events. Formally, these can all be thought of as Bremsstrahlung, putting appropriate incoming 



and outgoing lines in (11) (or a generalization of that equation to account for time dependent masses). 

An individual production, scattering, or decay event produces gravitational waves at all frequencies 
below the inverse timescale of the event. Multiple events, involving the same or different sources, can 
enhance the effect on the one hand, but also can introduce interference which suppresses the effect. 
For example, in a gas of particles with an approximately spherically symmetric distribution of particle 
positions and velocities, the quadrupole vanishes to first approximation and the net GW spectrum is 
a subleading effect. Moreover, if one considers the contribution of a single particle subject to multiple 
scattering events without relevant momentum loss, its net GW emission at very long wavelengths is 
simply determined by the scattering angle between the first and last event, with no enhancement from 



the additional events. This follows from the positive and negative contributions in (11) for incoming and 



outgoing lines. We will take these suppression factors into account in our estimates below. 
3.2.1 GWs from Production 

In this section, we discuss conditions under which the production event itself contributes a competitive 
tensor signal. Let us first consider gravitational waves, and then we will also address scalar emission. To 
analyze this, we need to specify the functional form of m(</>). During the production itself, we assume a 
coupling of the form <f >2 X 2 ^ so that m(4>) = <f> ~ <\>t. We will consider two examples for the later evolution: 
(i) m{(j)) continues to depend linearly on <p ~ <f>t, and (ii) m(0) transitions to a constant at some time 



t c < after the production event. In the appendix {A.l we describe a mechanism by which such 

a transition may arise. In a production event, the homogeneous rolling scalar field loses energy into x 
particles, and also into scalar radiation (reducing <fi in the process). 

Case (i): m((j)) = (ft « (f>t 
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Case (i) is interesting, as it is in some sense simplest to consider the (fi 2 x 2 model without assuming 
a more complicated functional form for the x mass. We will consider the GW emission arising from the 
sudden appearance of the produced particles and the associated scalar radiation modes. 



The stress-energy tensor for the particles, obtained by varying ( 13 ) with respect to the metric takes 
the form 

T p Z t = J2* (3) (Z-Ut)) P ^f0(t) (14) 
r> Pn 



where 



o 9<f>(t,x) _ g<j>(t,x)x 

V = i P = , = (15) 

V 1 — x 2 y 1 — x 2 



The spatial components are of the form 



TlJ ~ ^S(x - x n (t))9(t) ~ V ^tb{x - x n (t))8(t) (16) 



The additional i-dependence in the denominator in ( 16 ) translates into an additional factor of ui in its 
Fourier transform relative to the case studied in |12j which had a constant mass for the particles. The 
resulting gravitational wave emission is of the order 

dE GW fMW (17) 



du \M P J \EJ 



(This scales like the result |12j (11), times the extra factor of (co/E) 2 just noted.) 

From this we obtain at frequencies of order H that pew ~ Hn x H 2 /Mp where n x is the number 
density of produced particles and hence at frequencies uj ~ H we have 

,2 Pgw , H 3 

h ~ ~ j - , r r> (18) 

ptotal EM 2 



where we defined / = En x / ptotal (with ptotal ~ H 2 Mp) to be the fraction of the total energy density 
that is contained in the produced particles. As we explain more in detail in the next subsection, there is 
only a factor of n x here because all of the production events are independent. Turning this around, we 
see that 

where in the last step we put in the condition that h be detectible. The second factor here is > 1, 
so the inflationary scale H would have to be larger than 10~ e Mp in order for these particle sources to 
produce a detectible tensor signal. But this would require a higher scale of inflation than gives a standard 
contribution ([!]) to tensor modes of order h ~ 10~ 6 . 

However, there is also stress-energy in the inflaton field (f>, which can source GWs. We can estimate 



the scalar radiation as in the appendix, focusing on the upper limit of the integral in (98). This yields 
an energy density of order 

P5<t> ~ ^xtf 3 (20) 
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(at coupling g ~ 1), with the 5(f) particles carrying typical energies of order Eg^ ~ y 0. (The latter 
follows because we choose the upper range of k in (101), since we get the largest GW signal and the 
largest energy denstiy from the largest energy 5(f) particles.) 

From these 5(f) particles, we obtain GW Bremsstrahlung which is a factor of EjjH 2 times dl8h: 



h 2 ~ 



HEs,j, 



(21) 



Using that , with now Ess 



1 / 4 (i?M P ) 1 / 2 we get: 



f 2 /3,i/ 6 
J 56 6 ' 



/2/3 £ l/6 



x e 



1/6 



1/3 



(22) 



Here fss is the fraction of the total energy density which is carried by the 5(f) particles. This fraction 
is less than /, the fraction carried by the x particles, since those are ramping up in mass throughout a 

Hubble time, ultimately reaching energy E x ~ 4>/H. At g ~ 1 we have from (20) that fss ~ Hf/\ 



This was used in the last step of (22). Since H/Mp < h, and since we require 



j < e < 1, we see from 



2 X 2 coupling does not work at the level of production by 



this that the signal is still too weak to be competitive in this example 

Although this example with a simple 
itself, we will find below in §3.2.3 that subsequent decays can produce a very competitive signal. Before 
considering decays, however, let us consider a second model for m((f>). 

Case (ii): m — > const 

As our next example, let us consider case (ii) in which the mass becomes constant after a time 
t c <C H , and does not interact further in a Hubble time. In this case, for gravitational waves of 
frequency lo ~ H, we can work directly with the results (11) from the time- independent analysis of 
This gives 



h 2 ~ 



/ 



,EH 



(23) 



As above, we are interested in comparing the amplitude of gravitational waves produced by Bremsstrahlung 
to those produced by ordinary inflation Q, and so let us rewrite this as 



H 

~Mp 



h 



fhM P \ 



This means that for a given h the value of H/Mp needed to produce it is a factor of 

hM P 



(24) 



(25) 



smaller than the standard value. Since E < Mp and / < 1 we conclude it is only possible to reduce the 
Hubble parameter by an amount less than h which at the limit of detectability is roug hly 10~ 6 . Thus 
we have six orders of magnitude of potential gain. To obtain the full six orders requires having particles 



s 



of energy close to Mp having energy density comparable to that in the inflaton, but simply to obtain a 
competitive signal requires 

Eh 

w P >- j < 26 > 

The energy in our particles is of order their mass E ~ M. The asymptotic value of their mass depends 
on t c : 

M ~ g^t c , (27) 

or equivalently 

E 2 ~ M 2 ~ eg 2 Mp(Ht c ) 2 (28) 



Next we will study the consistency of the requirement ( 26 ) with the constraint imposed by the need 
to limit the scalar power emitted from the production event. The extent of the coupling to the scalar is 
determined by d ( / ) m(4>). Once the mass becomes constant, there is no longer any coupling to the scalar, 
but during the period between t = and t = t c in which the mass grows linearly with eft (and with 
t), there is a constant coupling to the scalar and we wish to estimate the scalar power from this whole 
process. 

For frequencies oo <C l/t c , there is destructive interference between the production event at t = and 
the event at t = t c when the coupling jumps to zero. This cancelation is exact at cv = 0, but for nonzero 
uj there is a residual contribution that arises from expanding a factor of e tujt that arises in the Fourier 
expansion of the radiation. In the radiated power, this introduces a suppression factor of (cot c ) 2 relative 



to the case analyzed in (107) where there was no time t c at which the scalars decouple. Because of the 
redshifting of modes within the Hubble patch before they freeze out, we get the largest contribution by 
taking modes of u ~ H and paying this (Ht c ) 2 suppression price. The scalar fluctuations are of order: 

C 2 ~^|(ift c ) 2 . (29) 
To assess the viability of the scenario we compute the scalar to tensor ratio, 

C 2 9 2 fH 2 M 2 1 

E iHtc) X jEH~ 7 2 - (30) 

Thus this scenario is viable as long as e ~ 10 _1 . 

Since the new source for £ fluctuations dominates, we need to ensure that the resulting scalar fluctua- 
tion satisfies the current bound on non-Gaussianities, at the level of 10 -3 . The amount of non-Gaussianity 
scales as l/y/N$, with representing the number of 5(f) fluctuations contained in an Hubble patch, and 
therefore we have the constraint > 10 6 . The constraint on the power spectrum gives an upper bound 
to the value of iV^, and we need to check that there is an open window. We have 

h 2 2 H 4 9 2 f 2 (Ht c ) 2 , , 



after using (p 2 ~ eH 2 Mj, and foty and (1281) to substitute for (H/M P ). We see that for extreme values 
Ht c ~ 1, g ~ 1, / ~ 1CT 1 and h ~ 10 -6 , can be as large as 10 10 , which gives quite a large 
window where the non-Gaussianity of the scalar fluctuations is compatible with observations, though 
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this allowed window shrinks if we move ourselves away from the most extreme region of parameter space. 
Interesting there is also a non negligible region of parameter space where the resulting non-Gaussian 
signal is detectable, but not yet ruled out. If all these constraints are satisfied, we can obtain a dominant 
tensor contribution from our sources, with H/Mp low enough to suppress the usual contribution ([!]). 

Finally, let us check that the scattering rate T is indeed < H, since otherwise we would need to take 
into account additional interactions, and that t c < H , since otherwise we would revert to case (i) above 
and pay an additional price in uj/E factors. The latter condition is (setting g ~ 1 for simplicity) 



Regarding the scattering rate T, we have 



an„v 



E < — 
H 



(32) 



£ 3 



The last factor v, the particle velocity, is p/yp 2 + m{t) 2 ~ p/m(t), where p 



the created particles. Now putting in (32) we see that 



H 2 



(33) 

is the momentum of 
(34) 



Since the latter factor is <C 1 in general in our model, this is consistent with a sufficiently slow decay 
rate, T < H. 



3.2.2 Decay of massive particles into massless ones 

In this section, we consider gravitational waves produced during decays of massive particles x present 
during inflation. As we have discussed, such particles may be produced via a coupling to the rolling 
inflaton such as g 2 4> 2 x 2 ] i n some circumstances [JjJ, this process repeats periodically during inflation. 
Here we assume that x couples to other light degrees of freedom Y, such that they can decay within 
a Hubble timescale. Similarly to what happens for the case of electromagnetic radiation in /3-decay, 
gravitational Bremsstrahlung radiation is produced not only in the case of scattering of particles, but 
also in a decay process. 

As we described above in {3.1, the amount of gravitational radiation per unit solid angle per unit 
frequency is given by [12J 

P D t ,D 



where 



F[{P N }, {P M }, k] = J2 P Tp .A (Pn ■ Pm? - \mW M ] ■ (36) 



N,M 



P N -kP M -k Ly 1V1 ' 2 



Here Di labels the collisions, x^ D . the location in space time where each collision/decay occurs, and {Pdj} 
is the set of momenta involved is each collision, rj is — 1 for incoming particles and 1 for outgoing ones. 
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For a single decay of one massive particle with mass M into two massless ones the above formula 
gives: 

/ rW, \ 1 M 2 

(37) 



/ dE \ 1 M 2 
\dujdn) ~ At^mJ 



\2 : r 



Here (E/M p i) is the effective dimensionless gravitational coupling squared, analogous to a in electro- 
magnetic Bremsstrahlung. 

For an arbitrary number of decays, the above formula becomes very complicated: events at different 
locations interfere if their relative distance is smaller than the frequency of the emitted wave, and the 



factor (36) involves momenta that belong to different events. In the present work, our main purpose 
is not to perform a precise calculation of the emitted gravitational radiation in one particular scenario. 
Most humbly we simply wish to estimate the order of magnitude of the effect in reasonable examples, in 
order to determine whether or not there is a direct relation between inflationary tensor modes and the 
height of the inflaton potential. 

Tensor modes at wavelengths much shorter than H~ l are highly suppressed by redshift (|6]). In 
each scattering, decay, or production event the energy produced in gravitational waves is frequency- 
independent for wavelengths longer than the timescale of the process. For these reasons, the leading 
contribution comes from those gravitational waves that are produced directly with Hubble frequency. 
Focusing on these wavelengths allows us to neglect the modulation due to the different location of the 
events, and consider them all at the same point in space. 

We must take into account possible destructive interference among the various decays, which result 



from the sum over momenta at different events in (36). The energy emitted in gravitational radiation goes 



as the square of the stress-energy tensor. So if there are N pait particles decaying in an Hubble patch, then 
naively the amount of gravitational radiation should go as N 2 &rt . Clearly this is an overestimate. If these 
particles are randomly distributed in an Hubble patch, their decay products will tend to be spherically 
distributed in the limit in which iV part is very large, and this will lead to a suppression of the emission 
of gravitational radiation. We can determine the net effect of the increased number of particles in the 
following way. The quadrupole-squared of the random distribution of particles will have a typical size 
proportional to -/V part instead of to N 2 t . This is very much like the variance of N independent random 
variables, which goes as N and not as N 2 , and also much like the typical realization of a random walk in 
quadrupole space, where we sum randomly all the quadrupoles associated to each event. Another way to 



think about this is to notice that the momentum sum in (36) assuming the two set of momenta belong 
to two different decays is not zero. However, it becomes zero if one averages over the outgoing direction 
of the momenta of the second decay. Now, in the case in which there are many decays happening at 



the same point, then the sum in (36) where the second momenta are taken from the various collisions 
effectively corresponds to averaging over the outgoing direction of the momenta of the second decay. In 
this case therefore the two sums over decays collapse to one, leaving us with a single factor of -/V par t. 

Altogether, if there are A^ part particles decaying within an Hubble patch, with A^ part large, the amount 
of produced gravitational radiation goes as 

( dE \ 1 M 2 

u^hi^Mj x (38) 



n 



This leads to the following amplitude in gravitational waves 



, 2 _ , 1 dp gw 



ptotoi din a; 



a;~_H" y total 



where n part is the number density of the decaying particles before they decay and where we have neglected 
numerical factors. If we call / the fraction of the energy density carried by the particles prior to their 
decay, we can re-write the above expression in a more useful way 

Let us see how big this number can be. 

Since the decaying particles do not redshift as approximately a cosmological constant, we need to 
have / < e, with e being the slow roll parameter, a number much smaller than one (but not necessarily 
tiny, let us say not larger than 10 _1 ). Notice that for standard slow- roll inflation, saturating this bound 
means that a fraction of order one of all the kinetic energy of the inflaton is dissipated in the creation 
of particles. In more general models this does not need to be the case, but this is an interesting regime 
to consider since it provides a mechanism to dissipate excessive kinetic energy on a steep potential, as in 
trapped inflation [9]. 

The mass M of the decaying particle can be bound to be at most of order M p i. We therefore can 
write the above formula as 

h 2 ~ e < e . (41 

M P M P e ~ M P K ' 

Alternatively, we can express the necessary value of H/Mp in order to have detectable signal: 

»m) , (42) 



JL h( hMp \ 

M P ~ \fM J 



which implies that the value of H/Mp can be reduced from the standard case by a factor of order 
hMp I \fM). This is quite a big improvement with respect to the ordinary case H/Mp ~ h. For example 
if we take / of order 10 -1 this means that we could detect gravitational waves for values of H/Mp as low 
as 10 -11 . This represents five orders of magnitude improvement with respect to the ordinary case (and 
ten orders of magnitude if you count in terms of the more physical parameter H 2 ). 

With a large window of opportunity for competitive GWs, let us consider more specific examples, 
relaxing some of the assumptions just made. For example, rather than taking M ~ Mp, we can consider 
the mass which is built up after the specific production mechanism in the appendix, arising from a 
coupling of the form g 2 <ft 2 x 2 , an d check how massive the x particles become before decaying. If we call 
At the time since the particle was created we can write: 



We can take (HAt) ~ 1 and get: 



M ~ g<t>At ~ g^Mp(HAt). (43) 
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If we consider / ~ e we obtain that H/Mp ~ 1CT 12 / ge 3 ^ 2 . If e ~ 1CT 1 and g ~ 1, we obtain H/Mp ~ 
3 x 1CP 10 , so the usual mechanism for tensor modes is suppressed by a large factor relative to the new 
sources in this example. 

We can be a bit more general and relax the (HAt) ~ 1 assumption. We can assume that there are 
many production events in a Hubble time, at a rate dNhits/dt and that in each of these production events 
n c is the number density of created x particles and that the decay rate of the particles is T^. The energy 
density in the x particles is given by: 



Px 



f dA t nc ^ e -(3^)A fM ^ n J^ptMl (45) 



with M ~ gy/eMp(Ht) and t ~ mm[l/H, l/T^]. We have also replaced the discrete sum over production 
events by an integral. We conclude that the typical mass of the x particles is determined by the shortest 
between a Hubble time and the lifetime because even if the lifetime is much longer than Hubble the 
abundance of very old particles dilutes exponentially. Only production events in the last t contribute 
particles at any given moment because particles from previous events have either decayed or diluted. 

We need to demand that the energy in the x particles and its decay products (Y) be a small fraction 
of the vacuum energy driving the expansion. We must take into consideration that the energy in the 
decay products of x does not have time to redshift during a Hubble time. This is especially relevant 
when Trf >> H. We have 

p x + p Y ~ fH 2 M P ~ n c ^Mi (46) 

which is just an energy M for each of the particles created in a Hubble time wether they have decayed 
or not. 

The amplitude of the tensor modes can be estimated by adding the contributions from all the x 
decays in a Hubble time. For simplicity we will approximate the density of x particles as constant during 



a Hubble time, given by p x in equation (45) divided by M. We then have: 



Mj, X \ 1, if T d >H, 140 

where the factor T^/H accounts for all the decays in a Hubble time. Thus for a fixed / the tensor 
amplitude is maximized when 1^ ~ H. For « H there are too few decays in a Hubble time to 
produce a lot of gravity waves and if >> H the mass of the x particles at the time of their decay is 
not large. We can turn this into an estimate for the Hubble scale for a given h: 

(48) 

This is the same result we got when assuming H At ~ 1 but enhanced by a factor H/Td or T^/H 
depending respectively on wether H is larger or smaller than T^. 

As before we need to make sure that the scalar power is below or at the same level as observed. The 
estimate is very similar to the example of gravity waves created at production which we analyzed before, 
in particular the case when the scalar coupling turned off. Again we could imagine than now the scalar 
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coupling turns off prior to the decay, at a time t c . If this is so the mass of a x particle does not continue 
to grow after t c and thus t c provides an upper limit for t, t — > min[i, t c \. As before the scalar power is 
suppressed by a factor (Ht c ) 2 . For the purpose of this suppression, the case when the scalar coupling 
never turns off corresponds to t c ~ t as the decay of the x particles with no scalar charge effectively 
acts in the same way as the shutting off of the scalar coupling and the Hubble time also provides a bound 
to the possible level of cancelations. Thus we always have t c ~ t. 
We can now compute the ratio between scalar and tensor power, 

ttt ^ x (Ht c ) 2 x — ~ -5, 49 

h 2 (M/Mp) 2 K ! 4> 2 e 2 V ' 

where the first term accounts for the ratio of couplings, the second for the suppression of Bremsstrahlung 
in the scalar case and the third comes from the conversion between <p and £ fluctuations and the normal- 
ization of the gravity wave energy density. Thus the scenario is viable as long as e ~ 10 _1 . Following the 



same steps that led to eq. (31), it is straightforward to check that the scalar power spectrum can easily 
satisfy the constraint on non-Gaussianity. 

Finally we need to check wether it is possible to neglect the annihilation of x particles into 5(p 
particles. The rate for this reaction goes like F ~ anv ~ aN part H 3 v (where v is their velocity and a their 
annihilation cross section). The cross section for 2 — > 2 scattering from the interaction term g 2 (j) 2 x 2 goes 
like 

a (8tt) 2 E 2 ~ (8tt) 2 lDUJ 

where in the last step we used that the energy E x carried by the massive x particles is > \J~(j) (it goes 

like M ^> \J~cf) at late times). Multiplying this through by N pSLVt H 3 v, and including a factor Nhu s for the 
number of production events during a Hubble time, we get 

r s H {■$?""•) (51) 

This is less than H, and hence completely negligible, as long as Nhn s < (87r) 2 /g |3 , a condition which is 
easy to satisfy. 



3.2.3 Creation of string pairs and their decay into rings of particles 

Next, let us consider a microscopic example leading to decay-induced Bremsstrahlung gravitational radi- 
ation. In this example, the rolling inflaton first produces pairs of strings, which decay into smaller string 
loops and then into particles. Consider a pair of long strings of length L produced via a time-dependent 
string tension, T ~ Tt. In the appendix below, we discuss the production of pairs of these strings, 
finding that the two members of the pair are created close to each other and with opposite orientation. 
These can quickly (on a timescale of T -1 / 3 <C H^ 1 ) fall apart into LT 1 / 3 smaller string loops J^] (Here we 
conservatively assume there are no small couplings in the system that suppress the interconnection of the 
strings; this is a feature of the "tensionless string theories" which arise on branes in the relevant string 

5 We thank J. Polchinski for this point. 
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constructions.) The total angular momentum of the original pair of long strings is zero, and the small 
loops into which it breaks do not have any preferred orientation or direction of spin; we will therefore 
treat these as random in our estimates. These smaller loops can then decay, in particular into scalar 
modes (and other light particles such as gravitons), with a random distribution of momentum directions. 
There is no quantum number protecting them from decay, and no small coupling which suppresses their 
decay rate. 

The original pair of long strings is not spherically symmetric, and supports a quadrupole. However, 
given the breakup into smaller loops and decay into light particles, for the purpose of long-wavelength 
emission we may describe the system as an instantaneous production of the final decay products. The 
final decay products, for example S(f> perturbations, do not have time-dependent masses, and so the 
analysis of p2] goes through unmodified. 



Including this, in the same way as in £3.2.2 we obtain 



dE T 2 / 3 



doodU M p 



2 ^loops-brings (52) 



where N r i ngs is the number of produced pairs of long strings (which then decay into rings of particles) 
Hence, using N loops ~ f 1 / 3 /H, 



£2/3 £1/3 £ 
PGw\u~H ~ H X H 3 X -^p—^ X brings ~ H 4 N rings (53) 



leading to 



h 2 „ P°K\^ H „ ( JL) 2 N ( T I (54) 
PTot ' \Mp) rmgs \HM 2 1 



In order to assess the strength of this effect, we need a model which determines T. In the scenario 
[SJ, T ~ r]Mp(p with rj < 1 a coupling. Since (fi/HMp ~ ^/e <C 1, the GW emission from a single ring of 
loops is subdominant to the contribution from tensor quantum fluctuations, by a factor of rj^/e <C 1. To 
get a competitive or better signal, we require 

Nrings > 

(55) 

We can bound N r i ngs above by imposing that the total energy NioopsNringsf 1 / 3 !! 3 be less than eH 2 Mp 
with N r i n g S > HMp/T so as to produce a viable signal. This leads altogether to the condition T 1 / 3 > H/ e, 
which is easy to satisfy. 

Finally, given the order one coupling of our strings to scalar perturbations 5cp, we must check if that 
would produce a contribution to the scalar power spectrum which is too large. We expect that the energy 
goes into scalars with wavelength of order T _1//3 to first approximation. In particular, decay into a large 
number of low-frequency modes is suppressed. The system is weakly coupled for energies much less than 
the string tension, which at time T~ 1 / 3 is of order T 2 / 3 . Emission into many 5<p particles of low frequency 
n is suppressed by a factor {oj/T l / 3 ) n an ". 

Given that the strings decay preferentially into 5<j) particles of frequency ~ T 1//3 , the resulting contri- 
bution to their energy density at freezeout is Aps^ ~ eH 2 Mp(H/T 1 / 3 ) 4: (where the last factor accounts 
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for the redshifting before freezeout). The ratio of this to the usual source of scalar perturbations (for 
which p 5(j) ~ H 4 ) is eH 2 Mp/f 4 / 3 . Evaluating this for T 1 / 3 ~ £H/e (with £ > 1), we obtain a ratio 
Apsfi/psfi ~ £,~ 4 e 5 Mp/H 2 . For e ~ 10~ 2 , this is easily subdominant for modest values of £. 

Finally, let us consider the possibility of Bremsstrahlung radiation into scalar modes. This is absent 
here: given the decay of the small loops into light particles such as 5cj) perturbations, the decay products 
do not couple to <j) and hence there is no appreciable scalar Bremsstrahlung. 



3.2.4 Bremsstrahlung from Scattering 

The final case we will study is the Bremsstrahlung radiation produced is when particles accelerate due 
to collisions. In a process where an energy E is transferred, Bremsstrahlung at low frequencies (u <C E) 
emits an equal amount of energy per unit frequency 

dE 9 ( E V 

^"UpJ • (56) 

for a single scattering event. However, in a gas of particles there will be multiple interactions. At 
sufficiently low frequency, the emission might be suppressed if particles interact before they can emit a 
graviton. We will parameterize this as: 



dE g ( E \ 2 f 1 , ^ w >7 2 T 



int 



du V Mp J I I ~r% — ) otherwise . 



(57) 



Here Yi n t is the rate of interaction in the frame in which our gas of particles overall has no net momentum 
and 7 is the boost factors of the scattering particles. 

These suppression factors [E] can be intuitively understood as the following]^] Consider a process 
where a particle interacts twice, hitting two targets which are a distance £ apart (in our cosmological 
frame). The rate of interaction is T ~ vn/£, where vn = v cos 9 is the component of the particle's velocity 
which is along the direction from the first to the second scattering event (in terms of an average scattering 



angle 6). The first interaction produces some GWs by Bremsstrahlung, and would give (56) if that were 
all that happened. In the second interaction, the particle emits by Bremsstrahlung a second graviton. 
This new wavefront interferes with the one emitted in the first interaction, effectively creating a higher 
frequency graviton. The frequency of the higher frequency graviton is determined by how far the first 
emitted wavefront gets before the particle scatters again. From the wavefront geometry we have 

uj- 1 =t/v\\ -l/c~lh 2 v\\ , (58) 

where in the second passage we have approximated 1 — v 2 /c 2 ~ 1 — v 2 /c 2 = I/7 2 which is valid for not 
too large deviation angles and v < c. We also have T = v cos8/£, so we can rewrite this in terms of T 
instead of i. This gives 

u ~ T 7 2 . (59) 



For massless particles, v = c and it is necessary to retain the 9 dependence in ( 58 ) , giving ui ~ Tcos(9) / ( 1 - 
cos{9)). 



3 Wc thank M. Pcskin for a useful discussion of this issue. 
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For the moment we take the factor of p as a free parameter, and will analyze the GW signal in various 
ranges of p. We will comment later on its possible values in particular cases which have been analyzed 
in the literature. It should also be stressed, as we will highlight next, that the threshold frequency 7 2 T 
at which the suppression starts, is only the first of a series of thresholds in which different physical 
mechanisms suppressing Bremsstrahlung become important and the value of p changes accordingly. We 
will see an example of this next. We take our initial threshold to be 7 2 r, which is the one that occurs 
when suppression is due to multiple scattering, as a guidance which is particular relevant for our setup. 
Generalization to different setups should be straightforward. 

Let us assume that there is a number density of particles n p each with typical energy E interacting 
with a rate Ti n t- Since during inflation we expect the density of gravitational waves to be stationary, we 
then have: 

^7^1Td7- (60) 

Even without assumption of stationarity, because of redshift, we are interested only in gravitons produced 
in about an Hubble time. This leads to the same factor of 1 / H above. By using n p = p p / E, and redshifting 
the gravitational waves down to the Hubble scale where they freeze out, we obtain: 

2 1 / dp gw \ Tint to dE g ( H 



h 2 — ^OL ^ /x l^ x _ x ^ x _ (61) 

Ptotai \dhiuj J ^ H H E dto \oj J 

where we defined / = p p / ptotai- 

We would like to understand if this is > 10~ 12 in a reasonable window of parameters, and whether 
this occurs in situations where the scale of inflation is too low to produce detectable tensor modes by the 
usual mechanism 0. 

Thermal equilibrium Since we are interested in gravitons produced in about an Hubble time, we 
need to have Ti n t to be at least as large as H. Since Ti n t is the numerator, let us start with a large 
Ti n t H. This implies we are in thermal equilibrium. We will include the possibility of a nontrivial 
species number and consider relativistic particles, giving 

p p ~ N*T A , E~ T , T int ~ N*a 2 T , (62) 

where a is the strength of the interaction. We can use this to solve for T: 



T ( f \ 1/4 / H \ 1/2 



M P \NJ \Mp) ' 
The amount of gravity waves produced in this case becomes 

2 H 3 T ( uj 



(63) 



int 



We will have to take into consideration the suppression of the Bremsstrahlung emission. Depending on 
the strength of the suppression, the value of the index p, it will be advantageous to consider waves emitted 
at either oj ~ H or to ~ 7 2 Ti n t- We consider each case in turn. 
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Index p > 3: In this case the suppression is strong enough that one is better off with GW with 
uj ~ 7 2 Ti n t. We then get: 

fc> ~'nS?- (65) 

Note that Ti n t is in the denominator so that we are better off having the smallest possible rate compatible 
with thermal equilibrium Ti n t ~ H. Also, in order to minimize the redshift, it will be convenient to take 
7 ~ 1. By using (63), we can rewrite (65) in terms of H/Mp to obtain: 



( — 

\M P 



h x 



/5 /6 




(66) 



The factor in parenthesis on the right represent how much we can reduce the value of H/Mp with 
respect to the standard value H/Mp ~ h. Reducing Ti n t/H to order one and by taking the extreme case 
~ / ~ 1 and h ~ 1CP 6 , we see that there are at most two orders of magnitude to be gained and 
the window closes as the interaction rate increases above H. More realistically / should be smaller than 
10 _1 , decreasing the possible improvement to about one order of magnitude. 

Index 1 < p < 3: In this case the suppression is sufficiently weak that one is better off with GW 
with uj ~ H. We have 



h 2 ~fx 



yi-P J1VT 
1 int n 1 

-f 2 PMp 



(67) 



Note that here with 1 < p < 3, 7 and Ti n t are in the denominator, so one better off with the lowest 
possible Ti n t ~ H consistent with thermal equilibrium and with 7 ~ 1. We obtain: 



(H_\ ( h l ^N l J G 
\Mp)~ h \ /5/6 




2(p-l)/3 



v 4p/3 



(68) 



which, after saturating the limit Ti n t ~ H and minimizing by taking 7 ~ 1, reduces to what we obtained 
the former case p > 3, leading to at most two orders of magnitude in possible gain, even in the extreme 
case. 



Index < p < 1 : In this case the suppression is not strong and again one is better off with GW 
with uj ~ H, but now Ti n t is in the numerator so we are better off with the largest possible Ti n t- Note 
that because 



h 1 



1 j 2 PM p 



(69) 



in the limit p — > there is a lot to be gained as the result becomes independent of H. However, we should 
consider the constraints coming from thermal equilibrium, as increasing Ti n t tends to increase T and this 



in turns tends to increase either / or H, both of which should be small. In fact, by using eq. (62) for 



T^t and then (63) for T, we obtain the improvement for generic p < 1: 




2+p 



2(l-p) 



(70) 
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This expression is minimized for / ~ a ~ 7 ~ 1. The maximum possible improvement would be given 
by P — > 0, leading to 

H , ( h \ 

which is about 6-7 orders of magnitude for reasonable number of species. However, p — > would mean 
that destructive interference is absent, which is not plausible. Determining the value of p for a given 



physical system (if the power law parameterization (57) applies) is quite complicated. 



In the world of collider physics, detailed studies of the Bremsstrahlung suppression have been per- 
formed in the case of a relativistic particle impacting on a fixed target, as this is the case of relevance in 
particle detectors and in fixed target experiments. Here we follow the recent treatment of |19j . where it is 
found that a value of p = 1/2 can be obtained in the case of ultrarelativistic scattering for frequencies in 
the interval 7 2 Ti nt > oj > ^int/l 2 - For lower frequencies oj < ^intjl 1 '', the suppression factor p becomes 
equal to p = 2 Q Notice that we have a non- vanishing interval of frequencies where the suppression is 



controlled by p = 1/2 only in the ultrarelativistic limit 7 > 1, so we cannot take 7 ~ 1 in (70), which 
would increase the possible gain. In this case (p=l/2), we have 



H ( ^5^/5 



M P \ fn/io N V w a 4/5 



(72) 



If we impose the constraint r > H > r/7 2 , then the gravitational modes oj ~ H which we are considering 
have frequency large enough to avoid the suppression p = 2, which kicks in at oj ~ r/7 2 , while still 
avoiding any suppression from redshifting. This contraint implies 

7 4/5 > ^3/10 /1/10a 4/5 (^V /5 >1, (73) 



By plugging the above expression into (72) and solving again for H/Mp, we find 

(74) 




which implies that even in this case there are at most two orders of magnitude to be gained. Plugging 
back on the constraint on 7, we obtain 

,>(^T>1, (75) 



h 

which can be realized achieving high boost 7 > 10 2 with reasonable number of species. 

'This is due to the fact that the small-angle approximation for the particle trajectory is no longer valid for 
such low frequency gravitons. We should point out that the minimum frequency at which we find Bremsstrahlung 
emission to be suppressed with a p — 1/2, the so-called LPM regime, is different from the value quoted in [T5] . 
But our crossover scale, determined by the geometry of the scattering events described above, agrees with the one 
quoted in the rest of the literature [20] . 



19 



It is not hard to imagine setups in which during the inflationary epoch there might persist the same 
conditions as in colliders, at least in principle. In particular, we need particles to interact in a highly 
boosted regime so that the scattering angle is small, since in the opposite regime the work [19] argues 
that the suppression factor turns into p = 2. As a proof of principle, let us take for example the case in 
which the inflaton produces two kind of species of particles due to its motion: one heavy non-relativistic 
scalar particle H and one light relativistic scalar particle L. If the only interaction is a quartic vertex 
HHLL, then we see that these light particle will thermalize by scattering with the bath of H particles, 
each interaction happening in the same kinematic regime as in cases studied for the particle colliders. LL 
elastic scattering is suppressed as it happens only due to the mediation of a loop of H particles, which 
can be thought to be very heavy with respect to the energy of the L particles. The Bremsstrahlung 
radiation produced in the LH interactions would then be expected to have a suppression factor p = 1/2, 
based on the results of [19] . 



Particles with a few additional interactions: So far in this section we have assumed that the 
particles were in thermal equilibrium, though sometimes we were led to take Ti n t ~ H. In this regime, 
it is not clear if the particles should obey a thermal distribution, and therefore, in order to explore the 
range of possibilities, it is interesting to consider the separate regime where the particle interact with a 
rate of order Hubble, but are not in thermal equilibrium. 

Consider the case when after production the particles interact a few more times in a Hubble time. 
This will allow us to solve for the energy using that: 

a 2 

H ~ T int = n p ((Tv} ~ n p —^, (76) 

with a strength of the interaction. This together with n p E ~ fptotal ~ fH 2 Mp results in an expression 
for the typical energy: 

which leads to 

h ~—^\M- p ) [iff) ■ (78) 
Depending on the value of p, we are better off with emitting gravity waves at frequencies of order 
7 2 r ~ -y 2 H {jp > 3) or H (p < 3). In either case we obtain: 

where q = 9/2 for p > 3 or q = 2>p/2 for p < 3. We can reduce H at most by three orders of magnitude. 

Away from equilibrium: At then end of this section, we are naturally led to consider the case in 
which Ti n t ^ H and thus we are not near thermal equilibrium. The particles we consider are a result 
of the decay of the inflaton so they go through at least one interaction, when they are created. The 
Bremsstrahlung from their creation process was analyzed in detail above in §3.2.1 and, as we saw, we 
obtained a larger window of GWs. 
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Scalar perturbations As in the previous examples, we must ensure that the scalar perturbations 
are consistent with current data and constraints. In the present case of production of gravitational waves 
through scattering, the connection to the scalar fluctuations is model dependent. We have analyzed the 
Bremsstrahlung radiation in this section assuming that the scattering particles do not have any coupling 
to the inflation, as occurs for appropriate mass functions. This reduces the question to whether the 
scalar modes from the production process are too great, in a case where the mass of the created particles 



becomes constant after production. As discussed above (29), we have 



C 2^HH {mc) 2 (80) 

e h/ 

Now t c is related to the mass M ultimately attained by the created particles: 

HM 1 M 

Ht c -. = — (81 

In the examples for which the scattering particles have 7> 1, the mass is much smaller than the energy 
E of the particles, so Ht c can be as small as we wish, suppressing £ 2 . 

In the examples with 7 ~ 1, the mass M must be of order the energy scale E of our sources. This 
leads to 

C 2 9 2 f(H\ gfHM 



e 3 / 2 \M P J y Cl e 2 M 2 



which needs to be < 10 10 . In the examples in thermal equilibrium above, using (|63l) we find 



With an order of magnitude gain in our GW signal, H/M p ~ 1(T 7 and this marginally fits. Similar 
comments apply to the case of particles with a few additional interactions. 



4 Discussion 

We have seen that a detection of tensor modes, even approximately scale-invariant tensor modes, from 
inflation does not automatically constitute a measurement of the inflationary potential energy. In this 
section, we discuss possibilities for distinguishing tensor modes from sources of the sort we have considered 
in this work from the standard source Q, in the event of a detection of primordial B modes. 

4.1 Power spectrum 

Our sources may be approximately scale invariant (e.g. in [HI [U]) with the inflaton coupling to new 
sectors of light degrees of freedom which are closely spaced along its trajectory. In that case, their power 
spectrum would be similar to that expected from ([!]). However, if the inflaton couples to light fields 
more sporadically, one could obtain a non-scale-invariant signal, which would distinguish the sources 
from 0. In the scale-invariant case, it would be intersesting to consider additional methods for breaking 
the degeneracy. 
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One approach is to consider non-Gaussianity. In the next subsection, we will discuss non-Gaussianity 
of the tensor spectrum. This was studied recently in |21j for vacuum fluctuations, where strong constraints 
on the shape were found. In the presence of a source sector of the kind we have here we expect a wider 
range of possible shapes for the three point function, though we will only estimate their magnitude here. 
It is worth noting that scalar non-Gaussianity arises in some of our examples: if the energy in produced 
sources is large (of order eH 2 Mp), and the mechanism [5j applies, then non-Gaussianity among the scalar 
modes is also predicted. 



4.2 A non-Gaussian window? 

In this section we have seen that there are several physical mechanism that could produce a signal in 
gravity waves that would be detectable and larger than the standard one due to the vacuum fluctuations. 
This means that we cannot derive the energy scale of inflation from observation of a scale invariant 
spectrum of gravity waves. It is worth pointing out that the scale invariance of gravity waves would still 
teach us about the time-translation invariance of the background when these modes were produced, a 
pristine signature of the quasi de-Sitter background which model-independently characterizes inflation 
[14] . It is however interesting to ask if it will be possible to distinguish between gravity waves as due 
to our mechanisms (in the scale invariant case) and the standard ones due to vacuum fluctuations. A 
possible distinction might arise from observation of the statistical properties of the gravity waves. Though 
it would be interesting to make a comprehensive study of this phenomenon using the Effective Field 
Theory of Inflation |14j , it is hard to imagine that vacuum fluctuations of gravitons during inflation can 
be very non-Gaussian. This is because the energy scale controlling the free Lagrangian for the gravitons 
is Mp, and this is a very high energy scale compared to the inflationary scale and to the usual canonical 
normalization of scalar fluctuations H l / 2 Mp/c s . Instead the new mechanisms of production that we have 
been describing in this section have nothing to do with vacuum fluctuations. Indeed, they are the result 
of interactions, and so are naively very non-Gaussian. However, there is another mechanism making to 
distribution Gaussian. This is the high number of gravitons being produced. Since each production event 
is independent, in the limit of high number of events the distribution becomes Gaussian, with a deviation 
of non-Gaussianity, parametrized by a dimensionless number that we call NG and that can be thought 
of as the skewness of the distribution rl that scales as the inverse square root of the number of gravitons 

— 1/2 

at a given frequency at horizon crossing -^g rav ^ ons * If is easy to estimate that the number of graviton is 
given by 

M 2 

^gravitons ~ H~ /5 gra vitons (w ~ H) ~ ~Jph ■ (84) 

Notice indeed that this is order one for vacuum fluctuations. If we define the gain factor <? ga in as a 
number in greater than one corresponding to the amount of decrease in H/Mp that we can have while 
still achieving a detectable signal with respect to the standard mechanism 

H h 

— (85) 

JWP Sgain 

we have 

(/gain ~ \/ -/Vg rav itons S'gain ~ JyQ • (^6) 



? In terms of the usual parameter /nl, we have NG ~ /nl^ 
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The amount of non-Gaussianity scales inversely to the gain. Observational constraints on NG scale as 

— 1/2 

the inverse square root of the number of modes that are signal dominated N , . Given the smallness 
of the signal in gravity waves, it is hard to imagine that non-Gaussianities will be tested below the level 
1CP 1 — 1CP 2 , implying that we would be able to distinguish these two scenarios only in the regime of 
modest gain, which is nevertheless a relevant fraction of the parameter space we have found with our 
examples. It is finally interesting to point out that there are also non-negligible regions of parameter space 
where the scalar fluctuations £ present some detectable non-Gaussian features. It would be interesting 
to study the details of the non-Gaussian distribution in the future. 
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A Appendix: Particle and String Production 

So far we have used the fact that a sector of sources X may be produced through time dependent motion 
of the inflaton. In this section, we will review and extend the computations of this effect. We assume 
couplings of the inflaton (ft to other degrees of freedom; this might be the generic situation and is necessary 
for reheating. For example, couplings of the form \g 2 (ft 2 x 2 endow the fields x with a time dependent 
mass 

m(tf = g 2 (ft{tf + ml (87) 

as (ft rolls. In string theory, it is as common to have time dependent string tensions, for example of the 
form 

T{tf = r] 2 Ml<f) 2 {t) + T 2 (88) 

for some dimensionless coupling rj that depends on the string coupling g s and the shape and size "moduli" 
of the extra dimensions. Sufficiently close to (ft = 0, we can approximate 

<f>(t) « vt (89) 

with v = eft the field velocity (of dimension 2). 

A simple way to see how such effects arise in string theory is to consider the realization of scalar fields 
from the relative motion of branes (as well as various dual descriptions of this). As two branes come 
together, strings stretching between them are particles in the worldvolume theory of the pair of branes 



which become light as in (87) (89). Similarly, membranes stretching between two branes constitute strings 
in the worldvolume theory. The tension T(t) of these strings is given by the membrane tension times the 
distance between the branes, which is proportional to (ft(t). (In both cases, if the branes miss each other 
a nonzero mass mo or tension Tq remains at the minimal distance between them.) In this section, we will 



rate of strings from (88) (89). 



first review the production of particles from ( 87 ) ( 89 ) and then generalize this to estimate the production 
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Figure 1: Light degrees of freedom can be particles (left figure) or higher dimensional defects such 
as strings (right figure). 

A.l More general mass functions 

For some of our examples in the main text, we wish to consider a mass m{(j)) for our produced source 
particles x which becomes approximately constant (independent of <fr) at some point after production. 
In this subsection, we describe a simple toy model with a few interactions depending on an additional 
heavy scalar (frn which yields this behavior classically upon integrating out (frn- This is closely related to 
a model described in [22] . 

The model has potential terms 

(M 2 H (cfr H - O ) 2 + 4> 2 H {4> - 4>i?) || = m(0, 4> H ?X 2 (90) 

where the index i refers to the production event which occurs at when the inflaton rolls through the point 
<j> = (f>% in field space. 

Once the particles are produced, there is a number density n x (which, as we will see below, is of order 
m 3 / 2 at the time of production) . This leads to energy density 

Px ~ m x n x ~ {X 2 )m(4>, 4>h) 2 - (91) 

Equivalently, 

m(<p,(pH) 



After the production event at <fr = (pi , at first (f> — (f>i is very small and the first term in ( 90 ) freezes (f>H 
at (fro. However, once <fr — (fri becomes as large as M#, the second term begins to dominate, and (frn adjusts 
to a more energetically favorable value (frn = 4>h* [</>l]) which depends on (fr^. Exactly as in the similar 
toy model of [22], the result of integrating out cfrn is that for <fr — fa S> Mh, the mass m((p, 4>h* [<^l]) 
becomes constant. 
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This type of model - and others more faithful to a concrete UV completion - may provide a rea- 
sonable starting point for studies of the radiative stability of the scenario where m((p) transitions to 
a constant. Incorporating supersymmetry helps suppress loop corrections, though this alone does not 
prevent Hubble scale mass corrections. Incorporating monodromy renders the corrections approximately 
periodic, restoring a discrete shift symmetry, but it is important to ensure that their amplitude is small 
enough to be viable. 

A. 2 Bogoliubov Coefficients: Particle case 

A particle with mass squared m(t) 2 = v 2 t 2 and frequency ui(t) = \/v 2 t 2 + k 2 has a WKB wavefunction 

^Jw-J—c-*/ 4 «(*')*' (93) 

which solves 

(-<9 2 - v 2 t 2 - h 2 )i>{t) = (94) 

to good approximation in the regime Co <C to 2 . 

We start with this pure positive frequency wavefunction in the far past and evolve it to the far 
future, where it picks up a negative frequency term proportional to the Bogoliubov coefficient f3 that 
determines the number of produced particles. In terms of creation and annihilation operators, this 
means a ear i y = aai ate + /3a/j at ■ So if we start in the vacuum, i.e. the state is the state killed by a ear i y , 
it is not killed by a\ a t e but instead is e^^ a ^ 2 ^ 2a \late vacuum > up to normalization. (We are describing 
it in the Heisenberg picture, for which states do not evolve.) Taking the expectation value of a\ at aiate 
reveals that |/3| 2 is the number of produced particles. 

In the above example, expanding co(t') ~ vt' + £17, the wavefunction at large t becomes, doing the 
integral in the exponent, 

ib(t) « -ie-^Tlfl 

Continuing t — > e~ l7T t, staying at large \t\ to preserve the WKB approximation [18J, pulls out the Bogoli- 
ubov coefficient 

Pk ~ e ' . 

Then the total number density of particles is J d?k\j3k\ 2 ~ v 3 ^ 2 . That is the dominant k is ~ y/v. 

This was done in the case with the mass going through zero. If instead oj 2 = v 2 t 2 + fj? + k 2 , we get 

|/3| 2 ~ e -*{k 2 W)h (95) 
So the number density of produced particles is model dependent, with v 3 / 2 being the upper bound. 



A. 2.1 Associated scalar emission 

In the bulk of this paper, we are interested in the gravitational wave emission from these produced 
sources. An important consideration is the level of scalar emission that accompanies their production 
and interactions. Here let us estimate this for the production event just reviewed. 
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A somewhat similar process was analyzed also in [9]. In the specific scenario worked out there, the 
scalar field dynamics is strongly affected by a very finely spaced set of points where different sectors of 
particles become light. In that case, because the scalar can lose energy to the many sectors of temporarily 
light fields, the Green's function for the scalar perturbation is not well approximated by its form in a free 
scalar field theory. For our purposes in the present work, we may consider a somewhat simpler regime in 
which these events are spread out enough that the scalar Green's function is the standard one. At the 
end of this section, we will verify that this is a self-consistent approximation. 

Let us calculate the scalar radiation emitted from the particle production event, treated as a version 
of Bremsstrahlung. This is similar to electromagnetic or gravitational Bremsstrahlung, but with the 

source in this case obtained from the Born-Infeld action — J drg^sj i 2 — x to be 

N x 

M*)=E / dr-S^^-x^r)) (96) 

n=l 3 7n 



where 7 n = i\ — ~x 2 l = ]P n l m n where dot is derivative with respect to r. From this, we obtain scalar 
radiation 



Ax '• d 3 k -m 2 n /p° n e- ik < x - x ^ 



^ = 'ZJw ~m >■> +c - c - (97) 

where we have summed over multiple x pairs produced in the event. 

We can simplify this by taking into account the fact that the particle slows down quickly; p° ~ m. 
Putting that in and considering a production event where N x particles are all produced at t = gives 

5^x)~g * e -*-*£ e *»». (98) 

3H k 71=1 

We have put the lower limit at H because below that we cannot use the flat spacetime analysis, and the 
upper limit at \J~(j) because the particle production process takes a time of order 1/\J~^> |^J Note that here 

k is the physical momentum, which we could call k/a{t) in terms of the comoving momentum k. 
Let us determine the expectation value of S(j)(x) 2 : 

Let us evaluate the term on the right inside the expectation value. Since particles are uncorrelated, this 
is proportional to 5 n ^ n > . This leads to 

e ik - £n Yl e~ ik '* n ') ~ ^(e i(fc - fc,) - f ") ~ N X H 3 / d 3 x e i{ - k ~ k '> Sn ~ N X H 3 (2vr) 3 5 {3) (k - k') , (100) 

n=l n'=l n=l 



9 In sec. 3.2.1 we consider the case in which x particles are coupled to <j> for a time shorter than an Hubble time, 



which, as we will see, leads to a suppression of the amount of 5(f> radiation produced at low frequencies. 
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where in the next to last step we have approximated the summation with the continuum limit and taken 
into account of the fact that the integral is limited to an order one Hubble patch. Plugging back into (99) 
and carrying one of the two k integrals through the delta function, we have 

(WWflXjP^st- (101) 

We are interested in the resulting curvature perturbation at horizon crossing, where £ ~ HS(j)/(f). Since 
fluctuations in 5(f) redshift as H/oj, with io being their frequency of emission, for each logarithmic frequency 
bin we obtain 

dP c H A 2 /H\ 3 , . 

—± ^9 N y — . (102 

dlogw x \uj J y ' 

Because of the redshift, the contribution to £ is dominated by 5(f) fluctuations emitted directly at u ~ H, 
leading to 

PC ~ ^ 9 2 N X < 1(T 10 , (103) 



where the last relation follows from the normalization. For large g 2 N x , we need a sufficiently small 
inflationary scale in order to match the normalization. This is an additional constraint, but one that can 
be satisfied for the examples given above. 

It is useful to obtain the same result as above in a way that is more similar to the way we obtain our 
result for gravitational waves. The energy density is 5(f) waves due to Bremmstrahlung is given by 

dps ndiSd) 2 ) dE d(5d> 2 ) n y dE 

H(p ~uj 2 Y ' ~n x uj— , => N Y ' ^— , (104) 

a log a; a log a; aw a log a; uj duo 

dE/du follows the same expression as for gravitational waves with (E/Mp) 2 replaced by g: 



' IL g\ (105) 



duj 
leading to 

d(5(p 2 ) g 2 n y g 2 H 3 N Y , . 

,\ ' ~ ~ K , 106 

a log uj oj uj 

which nicely agrees with what found above in (101). We can turn this into an estimate for ( 2 , 

a2 g 2 n x H 2 g 2 fH . . 

(107) 

and notice that this also implies 

Ns ~g 2 N x , (108) 

with Ns jX representing the number of respectively <f> and x particles in an Hubble patch. Finally, let us 
come back to the check we mentioned above, that we may use the free scalar Green's function in (97) 
in a useful range of parameters. To study this, we should compare the scalar perturbations 5(f) with the 
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distance A0 between particle production events. If the ratio 
do not typically lose energy to x production. By taking 5(f) ■ 



is less than 1, then the perturbations 
(5(j)(x) 2 ) and using (101), this ratio is 



HN h 



its 



N^ 2 H 2 



• N hlts < Pl /2 N h 



A(f> 



(109) 



where N^ts is the number of particle production events per Hubble time and where in the last passage we 



have used the constraint from (103) This ratio can be small consistently with a large Nhits, consistent 
with approximate scale invariance. (Because the events are discrete, there will be structure on small 
scales in the power spectrum, but for sufficiently finely spaced events these oscillations wash out.) 



A. 3 String Case 

It is interesting to generalize this to strings with a time dependent tension T(t) [7_]|^] In general this is 
not simply a sum over string oscillator states treated as particles [16], since the tension can vary rapidly 
enough that the string cannot causally adjust to maintain its oscillator configuration. We find that large 
pairs of loops are formed, but the relative oscillations between the strings and anti strings making up 
the pair are displaced from each other by a short distance, roughly of order T -1 / 3 . (This is required by 
causality.) They generically oscillate relative to each other, and may quickly decay into smaller loops 
of that size, though joining transitions forming longer loops also occurs at some level, depending on the 
density [17] . 



A. 3.1 Circular Loops 



Let us start by considering the simplest configuration (circular loops). The generalization of (94) for a 
circularly symmetric loop's wavefunction ^(r, t) is 



{-St + df - b l t l r l - k z )V(r, t) = 

where we took the time-dependent tension T(t) to behave as T(t) 2 = b 2 t 2 ; i.e. b = T. 
This can be derived from the Hamiltonian constraint in the string worldsheet theory 

dudr yf^T{t) {g^d a X M d p X N G M N) 



-dt 2 + dr 2 + r 2 d8 2 + dx 2 , 



Sstring ~ 

given in terms of the spacetime metric 



G M Ndx M dx N 



Varying (111) with respect to the metric produces the constraints 

= T a p : 



1 



T(t)(d a X M dpX n G MN - -g a ^X M d^X lv G M N) 



(110) 



(111) 



(112) 



(113) 
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One could use the upper bound of the integral in (101) to estimate S(f>/A(j) to ensure that even the most 
high-energy modes are not affected by particle production. This would lead to multiply the last term in ( 109 1 by 
a factor of (H 2 /^) 1 ^ 4 <C 1, which leads to an even milder constraint on Nhit s . 



li 



We thank J. Polchinski for very useful discussions. 



28 



We are interested in a simple circularly symmetric configuration, with 6 = a and r, x± and t being 
functions of r. For simplicity we are ignoring motion of the string in the r, 9 plane but this could be 
included. The T ar = constraint is solved by g UT = 0. The others are solved by 



i 2 + r 2 + 'x\ + r 2 



(114) 

Now writing this in terms of momenta, using Pt = = —T(t)i = idt (the last being the representation of 
the momentum operator in position space), and similarly for the other coordinates and momenta, yields 

{nob. 



In our problem (110), br is like v in (94) except there is also the d£\V term. However, if we work 
at r 3> t we have a regime where this term is subdominant in (110), basically because the derivatives 
with respect to r pull down factors of 1/r, which is smaller than \jt. That is, taking (93) with now 
uj 2 (t) = b 2 t 2 r 2 + k 2 , we get 

* « -oj 2 {t)^ = -{b 2 r 2 t 2 + k 2 )$> (115) 

but 

« -(bH 4 )W (116) 

For r ^> t, this is subdominant to the first, b 2 r 2 t 2 term in (3). It can also be subdominant to the 
k 2 term in (3) for the dominant k ~ ybr. In order for the large t expansion to be valid, we needed 
Co/oo 2 ~ l/(brt 2 ) <C 1. This seems to be consistent with the above analysis as long as r is large enough, 
r > fe -1 / 3 . 

According to this calculation, the time- varying tension T = bt can produce many pairs of large loops, 
because we get from this 

|/3| 2 „ e-ffcV&r (117) 

This was all done in the case that the tension goes through zero. One can also study similarly the case 
where it does not, giving (c.f. (95) with \x = T m i n r) 



ttk_ +T 2 . r ) 



(118) 



Although the loops can be large, the distance between them is mu ch sm aller, consistent with causality. 
Moreover, they move very slowly overall: the peak momentum from (117) is A;* ~ Mv ~ y/br. Using that 
M ~ b 2 ^r at the time t ~ of production, we find that the strings have a relative velocity 

(119) 



r i/3 J 



IS 



(120) 



&l/6 r l/2 

Therefore the time it takes to separate the two members of the pair by a distance T~ x / 3 

tseparation ~ b If/ 

This is much greater than the timescale T -1 / 3 ~ 6 -1 / 3 . 

This is important because the latter is the timescale for relative oscillations of the string in more 
general configurations. The circular configuration is not generic, and it is interesting to consider the 
problem more generally by studying the string path integral [7j. This yields a similar result, but with the 
two members of the pair of loops oscillating relative to each other at a relative distance 6" 1 / 3 ~ T -1 / 3 . The 
simplest string-theoretic examples which develop light strings (which are sometimes called "tensionless 
string theories") do not have a small coupling suppressing interconnection of the strings or decay into 
scalar modes. As a result, the pair production of large strings in these cases is quickly followed by the 
pair breaking up into a ring of smaller loops. This is the situation analyzed in §3.2.3 
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